Crossover from quantum to Boltzmann transport in graphene 
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We compare a fully quantum mechanical numerical calculation of the conductivity of graphene to 
the semiclassical Boltzmann theory. Considering a disorder potential that is smooth on the scale of 
the lattice spacing, we find quantitative agreement between the two approaches away from the Dirac 
point. At the Dirac point the two theories are incompatible at weak disorder, although they may 
be compatible for strong disorder. Our numerical calculations provide a quantitative description of 
the full crossover between the quantum and semiclassical graphene transport regimes. 
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Arguably, one of the most intriguing properties of 
graphene transport is the non-vanishing "minimum con- 
ductivity" at the Dirac point. The carrier density n in 
these single monatomic sheets of carbon can be continu- 
ously tuned from electron-like carriers for large positive 
gate bias to hole-like carriers for negative bias The 
physics close to zero carrier density (also called the intrin- 
sic or Dirac region), is now understood to be dominated 
by the inhomogeneous situation where the local poten- 
tial fluctuates around zero, breaking the landscape into 
puddles of electrons and holes Q. In the literature, two 
separate pictures have emerged to understand the value 
of this minimum conductivity. The first picture expands 
around the universal value for the minimum conductivity 
frnin = 4e 2 / (71-/1) for clean graphene 0] and argues that 
the presence of potential fluctuations smooth on the scale 
of the graphene lattice spacing increases the conductiv- 
ity through quantum interference effects 0, H, [f| ■ The 
second picture extrapolates to the Dirac region from the 
high-density limit, where the conductivity for charged 
impurity scattering can be calculated using the semiclas- 
sical Boltzmann theory. This approach has been further 
refined near the Dirac point by positing that the sys- 
tem acquires an effective carrier density n* calculated 
from the rms density fluctuations (associated with the 
electron-hole puddles) about the Dirac point caused by 
the same impurities that are responsible for the scattering 
of carriers at high density [|| . These two conceptually dif- 
ferent approaches lead to strikingly different predictions 
for the conductivity at the Dirac point: The Boltzmann 
transport theory for Coulomb disorder predicts that in- 
creasing disorder decreases the conductivity, whereas the 
weak antilocalization picture has a conductivity that in- 
creases with increasing disorder strength. Given their 
vastly different starting points, it is not surprising that 
the two approaches disagree. 

A direct comparison between the two approaches has 
not been possible, mainly because the published pre- 
dictions of the Boltzmann approach include screening 
of the Coulomb disorder potential, whereas the fully 



quantum-mechanical numerical calculations are for a 
non-interacting model using Gaussian disorder. Notwith- 
standing the fact that screening and Coulomb scattering 
play crucial roles in transport of real electrons through 
real graphene [§] , the important question of the compar- 
ison between quantum and Boltzmann theories has re- 
mained unanswered, even for the Gaussian disorder case. 
It is the goal of this work to provide such a comparison, 
thereby establishing the bridge between these two widely 
used complementary theoretical approaches to transport 
in graphene. 

In what follows we consider non-interacting electrons 
at zero temperature in a Gaussian correlated disorder po- 
tential that varies smoothly on the scale of the graphene 
lattice spacing. This situation is described by the effec- 
tive Hamiltonian H = a ■ p + U(r), where v is the Fermi 
velocity, p the (two-dimensional) momentum, and U(r) 
a random Gaussian potential with correlation function 
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where £ is the correlation length and Kq is a dimension- 
less parameter that parameterizes its magnitude. Typical 
experimental conditions correspond to Kq between 1 and 
3 [lfj, [ll ] . We numerically solve the full quantum prob- 
lem for a sample of finite size L>(, starting from intrin- 
sic graphene where quantum coherence effects dominate, 
to high doping where quantum effects are a small cor- 
rection to the conductivity fl2| . We have compared the 
numerical results to predictions of the Boltzmann theory, 
its self-consistent modification of Ref. d, and weak antilo- 
calization corrections, for a range of disorder strength Kq 
and carrier densities n. 

Our main conclusions, to be supported by the mate- 
rial below, are: (i) Away from the Dirac point, both the 
Boltzmann theory and the full quantum solution agree to 
leading order, a cx n 3 / 2 , with deviations only in terms 
of order n 1 / 2 and smaller. This validates the assump- 
tions of both theories, i.e., the Born approximation for 
the Boltzmann approach and use of a finite sample size 
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in the numerical approach, (ii) At the Dirac point, the 
quantum conductivity increases with increasing disorder 
strength K ; for K » 1 the increase is compatible with 
the self-consistent Boltzmann theory, (iii) As a function 
of carrier density, the quantum conductivity, but not the 
Boltzmann conductivity, shows a sharp reduction at the 
Dirac point, which is most pronounced for Kq ~ 1; the 
conductivity becomes proportional to 1/Kq away from 
the Dirac point; (iv) The numerical quantum results are 
consistent with da/dhiL — 4e 2 /irh for a > 4e 2 /h, ir- 
respective of the carrier density n, consistent with the 
weak antilocalization theory. Our numerical calculations 
provide a quantitative description of the full crossover be- 
tween the quantum and semiclassical transport regimes, 
for which no analytical theory is available. 

The Boltzmann conductivity corresponding to the 
model of Eq. Q] is calculated using the relation a = 
e 2 vv 2 7 r/h, where v = Akp/^tiVF is the density of states 
and the elastic relaxation time r is given by 
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where # q parameterizes the direction of q, so that 
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kp/ir. The leading term for 



with the carrier density n 
large density can also be obtained considering the clas- 
sical diffusion of a particle undergoing small-angle de- 
flections from the random potential U [HI . The weak 
antilocalization correction to the conductivity is [i[ 



Sa(L,£) = In (L/£).. 
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where £ is the transport mean free path. In the Boltz- 
mann theory, £ can be obtained from the relation <jb = 
2(e 2 /h)kp£. A self-consistent modification of the Boltz- 
mann theory was proposed in Ref. [g in order to describe 
transport near the Dirac point n — 0. For our Gaussian 
model of disorder, this modification involves replacing 
the carrier density by a "self-consistent" carrier density 
n* = ir^iep/hvp) 2 , where ef = {{e F + U) 2 ) Q. We 
then find n* — \n\ + Kq/2tt 2 ^ 2 , and the self-consistent 
prediction for the conductivity is given by Eq. ([2]) above 
with n replaced by the self-consistent density n*. 

In the numerical calculation we consider a graphenc 
strip of dimensions L x W with W, L ^> £ , connected to 
a highly doped graphene regions on both ends. Following 
the method described in Ref. [5|, we calculate the con- 
ductance G of the graphene strip. The conductivity a is 
then obtained using the relation 



R = l/G. 
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FIG. 1: (Color online) Resistance, R — l/G (left) and con- 
ductivity a obtained using Eq. @ (right), as a function of 
sample length L. The three curves shown are for W/£ = 200, 
Ko = 2 and 7m£ 2 = 0, 0.25, and 1 [from top to bottom (bot- 
tom to top) in left (right) panel]. Dashed lines in the right 
panel show da/ din L = 4e 2 /7r/i. The inset in the left panel 
shows the crossover to diffusive transport (i > Q. 



We verify that our results do not depend on the real-space 
discretization in the longitudinal direction, the cut-off of 
the transverse momentum (see Ref. [| for details), and 
the aspect ratio W/L. Extracting the conductivity using 
((31) is different from Ref. 0, where the conductivity was 
identified with LG/W. The advantage of Eq. (0]) is that 
it eliminates the effect of an additive resistance from a 
region of ballistic transport adjacent to the contacts to 
the source and drain reservoirs and, hence, gives accurate 
conductivities for smaller samples sizes than the identifi- 
cation of a and LG/W . Our procedure is illustrated in 
Fig. [I] where we show typical quantum numerical results 
for the resistance R = l/G and the conductivity <r(L) 
defined through Eq. ([4J. 

We restrict the analysis of our numerical data to sam- 
ples with length L > £. (In semiclassical transport, this 
is the regime where electron motion is diffusive.) Ac- 
cording to the Boltzmann theory, £ ~ £ 3 n/A"o- In the 
diffusive regime, the quantum conductivity a has a weak 
dependence on the sample length L because of weak an- 
tilocalization. In order to compare numerical data at 
different Kq or n, we use two different procedures. In 
Figs. [2] and [31 we compare conductivities at a reference 
sample length L — 50£, which is well inside the diffusive 
regime if Kq > 1. For this sample size and beyond, the 
L-dependence of the conductivity was consistent with the 
theoretical expectation da/d\nL = Ae 2 /irh of weak an- 
tilocalization theory for all carrier densities if Kq > 1. 
(Representative data for Kg — 2 shown in Fig. [TJ) At 
the Dirac point, this size-dependence of the conductivity 
was previously observed in Refs. H, 0, EH- 

Figure [2] shows the conductivity <r(L — 50£) versus 
the carrier density n for two values of Kq. For n£ 2 > 1 
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FIG. 2: (Color online) Conductivity a versus carrier density 
n. The left panel shows the low-density behavior near the 
Dirac point, the right panel shows the high-density behav- 
ior. Data points are from numerical simulation with Ko = I 
(squares) and Ko — 4 (diamonds) with dotted lines in the 
left panel as a guide to the eye. The dashed curves show the 
predictions of the Boltzmann theory, and the solid lines show 
the self-consistent Boltzmann result. 



the conductivity is well described by the asymptotic be- 
havior of Eq. @, the dominant correction term being 
proportional to n 1 / 2 . Replacing the carrier density n by 
the self-consistent carrier density n* (solid lines) further 
improves the agreement. For small densities, the quan- 
tum conductivity shows a sharp minimum at the Dirac 
point n = 0, which is most pronounced for small disorder 
strengths. Such a dip is not present in either the Boltz- 
mann theory or its self-consistent modification. Since in 
the quantum theory a increases with increasing Ko at the 
Dirac point but decreases with increasing Ko away from 
the Dirac point, the quantum a vs. n curves of different 
Ko cross somewhere in the region < 7m£ 2 < 1 for the 
parameter range we consider. This reversal in behavior 
was previously noted by Lewenkopf et al. in numerical 
simulations of a tight-binding model (lij . although the 
numerical data of Ref . [IH do not allow a conclusion to be 
made about large carrier density. The agreement at high 
carrier density between the quantum and the Boltzmann 
theory is an important new result of this work. 

In Fig. [31 we address the conductivity as a function of 
disorder, comparing the quantum and Boltzmann theo- 
ries. Motivated by prediction of Ref. d that the a vs. 
n curve exhibits a plateau of width ~ Kq/£, 2 near the 
Dirac point, we consider the conductivity at the Dirac 
point n — (left panel) as well as near the edge of 
the proposed plateau, at mi = Ko/(irt;) 2 (right panel) 
171 ] . The numerical calculations at the plateau edge are 
in good qualitative agreement with the self-consistent 
Boltzmann theory. At the Diract point, however, a is 
found to increase with Ko for the entire parameter range 
we consider, which differs from the prediction of the 
Boltzmann theory 18J and the self-consistent Boltzmann 
theory. The former predicts a = 8e 2 /Koh at the Dirac 




FIG. 3: (Color online) Conductivity a versus disorder 
strength at the Dirac point (left) and at carrier density 
irn = Ko/(n^) 2 , corresponding to the edge of the minimum 
conductivity plateau of Ref. 8 (right). Data points are from 
the numerical calculation for L = 50£ and the (solid) dashed 
curves represent the (self-consistent) Boltzmann theory. 



point, whereas the latter deviates from this prediction for 
Ko ~ 1, reaches a minimum at Ko ~ 9.71, and crosses 
over to the asymptotic dependence a ~ 2e 2 K^ 2 /irh for 
Ko ^> 10. At large Ko the numerical data follow the 
trend of the self-consistent theory, although we cannot 

1/2 

confirm the asymptotic dependence oc K from the pa- 
rameter range studied in our simulations. Upon reducing 
Ko below unity, the conductivity first decreases sharply, 
consistent with a renormalization of the mean free path 
£ for Ko < 1 @, Upon reducing Ko further, the 

Dirac point conductivity saturates at the ballistic value 
(7 = 4e 2 /nh. 

The system-size-dependent weak antilocalization cor- 
rection ([3]) is included in Figs. [2] and [3j which show a 
at the reference length L — 50£. In Fig. 0] we subtract 
the L-dependent logarithmic increase and show the Ko 
dependence of a' = ]imi / - too [a(L) — n^ 1 ln(L/£)] [20| . 
Subtracting weak antilocalization significantly improves 
the agreement with the self-consistent theory at large Ko- 
Unlike the conductivity at the reference length L = 50£, 
which saturates at 4e 2 /irh for small Ko, o' continues to 
decrease without bounds if Ko is lowered. 

The increase of a with Kq at the Dirac point for weak 
disorder is markedly different from the prediction of the 
Boltzmann theory. A key assumption of this theory and 
its self-consistent modification is that the graphenc elec- 
tron liquid can be mapped to an essentially homogeneous 
system with an effective carrier density n* equal to the 
rms of a fluctuating "local" density determined by the 
random potential U. This assumption becomes question- 
able at the Dirac point, where the electron liquid is bro- 
ken up in puddles of electron-like and hole-like regions. 
At weak disorder, Ko <C 1, quantum fluctuations spread 
the carriers over many puddles and the concept of a local 
carrier density becomes problematic. It is in this regime 
that the difference between the quantum and Boltzmann 
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FIG. 4: (Color online) Same as Fig. but for a' = 
]im L ^oo[<r{L)--K- 1 hx(L/0]- 



calculations are, as expected, most pronounced. 

The Gaussian random potential used here is the 
potential-of-choice for comparisons of analytical theories 
and numerical simulations. Yet, it differs in essential 
ways from the random potential in realistic graphcnc 
samples that do not have Gaussian statistics, since it is 
likely caused by charged impurities in the substrate with 
a typical distance d from the graphcnc sheet smaller than 
the spacing between impurities. Still, it may be possible 
to extract equivalent parameters Kq and £ from a real- 
istic random potential (see Ref. [Tlh . implying that the 
sharp dip in conductivity predicted in the quantum the- 
ory would occur in a window n ~ 5 x 10 10 cm~ 2 around 
the Dirac point 1211 . T his feature has not been observed 
in experiments jlO ,122, 2^, 24 \. Reasons why the dip has 
not been seen could be a suppression of quantum coher- 
ence by finite temperature effects or rippling [25| of the 
graphcnc sheet, or long-range fluctuations of the mean 
carrier density which effectively smear the feature near 
n = [H. 
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